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A  theoretical  model  for  prediction  of  burning  rates  of  composite 
(ammonium  perchlorate  oxidizer)  solid  propellants  as  a  function  of  pressure 
and  crossflow  velocity  has  been  developed.  Included  in  this  model  is  the 
capability  for  treatment  of  multimodal  oxidizer  particle  sizes  and  metalized 
formulations.  In  addition,  an  experimental  device  for  measuring  the  effects 
of  crossflow  velocity  on  propellant  burning  rate  has  been  developed  and  used  — 
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to  characterize  a  series  of  AP/HTPB  propellants  with  systematically  varied 
formulation  parameters.  Model  predictions  of  zero-crossflow  burning  rate 
versus  pressure  characteristics  have  been  found  to  be  in  excellent  agreement 
with  data,  while  the  agreement  between  erosive  burning  predictions  and  data 
is  in  general  good.^  The  experimental  data  indicate  that  the  dominant  factor 
influencing  the  seJifivLtivity  of  composite  propellant  burning  rate  to  crossflow 
is  the  base  (no-crossflow)  burning  rate  versus  pressure  characteristics  of  the 
propellant  (lower  base  burning  rate  leading  to  increased  crossflow  sensitivity) 
with  other  factors  having  at  most  a  second  order  effect  outside  their  influ¬ 
ence  on  base  burning  rate.  For  example,  three  formulations  with  widely  diff¬ 
erent  compositional  and  ingredient  particle  size  parameters  but  with  nearly 
identical  base  burning  characteristics  exhibited  very  similar  erosive  burning 
characteristics.  Finally,  the  model  has  been  used  to  examine  the  effects  of 
motor  scaling  on  erosive  burning:  erosive  burning  is  predicted  to  diminish 
with  increasing  motor  size,  in  agreement  with  experience. 
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PREDICTED  AND  MEASURED  EFFECTS  OF  PRESSURE 
AND  CROSSFLOW  VELOCITY  ON  COMPOSITE 
PROPELLANT  BURNING  RATE* 


Merrill  K.  King 
Atlantic  Research  Corporation 
Alexandria,  Virginia  22314 

ABSTRACT 

A  theoretical  .model  for  prediction  of  burning  rates  of  composite 
(ammonium  perchlorate  oxidizer)  solid  propellants  as  a  function  of  pressure 
and  crossflow  velocity  has  been  developed.  Included  in  this  model  is  the 
capability  for  treatment  of  multimodal  oxidizer  particle  sizes  and  metalized 
formulations.  In  addition,  an  experimental  device  for  measuring  the  effects  of 
crossflow  velocity  on  propellant  burning  rate  has  been  developed  and  used  to 
characterize  a  series  of  AP/HTPB  propellants  with  systematically  varied 
formulation  parameters.  Model  predictions  of  zero-crossflow  burning  rate 
versus  pressure  characteristics  have  been  found  to  be  in  excellent  agreement 
with  data,  while  the  agreement  between  erosive  burning  predictions  and  data  is 
in  general  good.  The  experimental  data  indicate  that  the  dominant  factor 
influencing  the  sensitivity  of  composite  propellant  burning  rate  to  crossflow 
is  the  base  (no-crossflow)  burning  rate  versus  pressure  characteristics  of  the 
propellant  (lower  base  burning  rate  leading  to  increased  crossflow  sen¬ 
sitivity)  with  other  factors  having  at  most  a  second  order  effect  outside  their 
influence  on  base  burning  rate.  For  example,  three  formulations  with  widely 
different  compositional  and  ingredient  particle  size  parameters  but  with 
nearly  identical  base  burning  characteristics  exhibited  very  similar  erosive 
burning  characteristics.  Finally,  the  model  has  been  used  to  examine  the 
effects  of  motor  scaling  on  erosive  burning:  erosive  burning  is  predicted  to 
diminish  with  increasing  motor  size,  in  agreement  with  experience. 

INTRODUCTION  AND  BACKGROUND 

Erosive  burning,  the  alteration  of  propellant  burning  rate  by  high 
velocity  product  flow  across  the  burning  surface,  has  become  an  increasingly 
important  phenomenon  with  the  advent  of  very  low  port-to-throat  area  ratio 
cylindrically  perforated  motors  and  nozzleless  motors.  The  motor  designer 
must  be  able  to  predict  this  burning  rate  modification  and  in  particular 
understand  the  effects  of  motor  scaling  and  length/diameter  ratio  on  erosive 
burning  in  order  to  properly  carry  out  his  function.  In  addition,  the 
propellant  chemist  needs  to  understand  the  effects  of  various  formulation 
parameters  on  the  sensitivity  of  a  propellant  to  crossflow  in  order  to  tailor 
propellants  to  desired  ballistic  characteristics. 


*Research  sponsored  by  the  Air  Force  Office  of  Scientific  Research  (AFSC), 
United  States  Air  Force,  under  Contract  F49620-78-C-0016.  The  United  States 
Government  is  authorized  to  reproduce  and  distribute  reprints  for  Governmental 
purposes  notwithstanding  any  copyright  notation  hereon. 
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Over  the  past  several  years,  this  investigator  has  been  conducting  an 
experimental  and  analytical  study  o£  the  erosive  burning  of  composite 
propellants  with  the  aims  of  determining  how  various  formulation  parameters 
influence  the  sensitivity  of  propellant  burning  rate  to  crossflow,  developing 
an  a  priori  burning  rate  model  (second  generation  model)  which  will  permit 
accurate  prediction  of  composite  propellant  burning  rate  as  a  function  of 
pressure  and  crossflow  velocity  given  only  the  formulation  composition  and 
solid  ingredient  particle  sizes,  and  developing  scaling  laws  to  permit 
extrapolation  of  erosive  burning  data  obtained  in  test  devices  and  small  motors 
to  larger  motor  situations.  Details  of  the  model  development,  with  the 
exception  of  the  treatment  of  effects  of  aluminum  additive  on  burning  rate,  are 
presented  in  References  1  and  2.  During  the  past  year,  the  model  has  been 
slightly  modified  to  more  accurately  treat  the  augmentation  of  transport 
properties  by ‘crossflow-induced  turbulence,  and  it  has  been  extended  to  treat 
aluminized  formulations  as  well  as  non-metalized  ones. 

Since  there  is  a  dearth  of  systematic  experimental  erosive  burning  data  in 
the  literature,  a  test  device  has  been  developed  and  used  to  characterize  the 
erosive  burning  behavior  of  a  series  of  AP/HTPB  composite  propellants  with 
systematically  varied  formulation  properties.  Zero-crossflow  predictions  and 
data  have  been  generated  for  thirteen  non-catalyzed  AP/HTPB  formulations  (four 
containing  aluminum)  while  crossflow  data  and  predictions  have  been  obtained 
for  ten  of  these  formulations  to  date.  A  brief  description  of  the  model  (with 
emphasis  on  its  extension  to  treat  metalized  propellants),  comparison  of 
theoretical  predictions  with  data,  and  definition  of  major  parameters  af¬ 
fecting  composite  propellant  crossflow  sensitivity  are  presented  in  the 
following  sections. 


MODEL  DESCRIPTION 

The  major  elements  involved  in  development  of  a  complete  a  priori  model 
for  prediction  of  composite  propellant  burning  rate  as  a  function  of  pressure 
and  crossflow  velocity  include  development  of  a  no-crossflow  composite 
propellant  model  (embodying  many  of  the  concepts  of  the  Beckstead-Derr- 
Price^)  model)  for  prediction  of  burning  rate  as  a  function  of  pressure  in  the 
absence  of  crossflow,  followed  by  incorporation  of  modifications  based  on  two 
postulated  mechanisms  for  augmentation  of  heat  feedback  from  flames  in  the  gas 
phase  above  the  propellant  surface  by  crossflow.  (See  Figure  1).  As  discussed 
in  Reference  2,  several  variants  of  a  basic  no-crossflow  model  were  developed 
for  unimoda 1-oxidizer  non-metalized  propellants,  with  one  variant  finally 
being  selected  for  extension  to  treatment  of  formulations  containing  multi¬ 
modal  oxidizer  and  metal  fuel. 

The  basic  model  centers  around  an  energy  balance  at  the  propellant 
surface.  In  this  balance,  the  product  of  the  propellant  burning  mass  flow  and 
the  heat  required  to  raise  the  ingredients  from  ambient  temperature  to  the 
surface  temperature  (related  to  the  burning  rate  by  an  Arrhenius  function)  plus 
the  heat  required  to  vaporize  that  fraction  of  the  ingredients  which  do  not 
exothermically  react  just  below  the  surface  is  equated  to  the  sum  of  the  heat 
release  rate  from  subsurface  reactions  and  the  rate  of  heat  feedback  from  the 
two  gas  flame  zones  depicted  in  Figure  1.  (For  metalized  formulations,  as 
discussed  later,  an  additional  heat  release  zone  associated  with  burning  of  the 
metal  in  the  gas  above  the  propellant  surface,  not  depicted  in  Figure  1,  is  also 
considered.)  Thus,  the  burning  rate  of  non-metalized  formulations  is 
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controlled  by  three  principal  heat  release  zones:  (1)  heat  release  in  a  thin 
subsurface  zone  just  adjacent  to  the  propellant  surface;  (2)  heat  release  in 
the  gas  piiase  above  the  propellant  from  ammonium  perchlorate  decomposition 
products  burning  as  a  nonpropellant;  and  (3)  heat  release  from  a  diffusion 
flame  between  the  AP  decomposition  (and  monopropellant  flame)  products  and 
fuel  vapor  released  by  binder  pyrolysis. 

The  subsurface  heat  release  is  calculated  by  an  iterative  process, 
coupled  with  the  remainder  of  the  model,  in  which  an  estimate  of  the  subsurface 
temperature  profile  is  made  and  substituted  into  an  Arrhenius  rate  expression 
representing  subsurface  heat  release  data  measured  by  Waesche  and  Wenograd^), 
which  is  then  integrated  from  the  surface  to  a  depth  where  the  temperature 
drops  below  the  melting  point  of  AP  to  obtain  the  total  subsurface  heat  release 
per  unit  mass  of.  propellant.  This  is  then  multiplied  by  the  burning  mass  flux 
to  yield  a  heat  release  rate.  This  procedure  differs  markedly  from  that  of  the 
BDP  model-  in  which  the  amount  of  subsurface  heat  release  per  unit  mass  of 
propellant  is  assumed  to  be  a  constant,  independent  of  such  parameters  as 
burning  rate,  and  is  included  with  the  binder  heat  of  vaporization.  Since  the 
subsurface  temperature  profile  steepens  rapidly  with  increasing  burning  rate, 
while  surface  temperature  increases  with  burning  rate,  the  procedure  used 
herein  results  in  the  subsurface  heat  release  per  unit  mass  of  propellant 
varying  with  the  burning  rate. 

For  the  gas  phase,  a  two-flame  approach  was  chosen  for  this  model,  the  two 
flames  being  an  AP  monopropellant  flame  and  a  columnar  diffusion  (Burke- 
Schumann^))  flame.  As  indicated  in  Figure  1,  three  distance  parameters 
(FH9Osin0,  L^p,  and  Lpy)  are  important  in  calculating  heat  feedback  from  these 
gas  flames  to  the  propellant  surface.  FH90  refers  to  the  distance  associated 
with  completion  of  mixing  of  90  percent  of  the  fuel  and  oxidizer  gas  products, 
while  Lpx  and  L^p  are  reaction  distances  (products  of  reaction  times  and  gas 
velocity  away  from  the  surface)  associated  with  the  binder  gas-oxidizer  gas 
flame  and  the  monopropellant  AP  gas  flame,  respectively.  As  discussed  below, 
flame  bending  associated  with  crossflow  is  postulated  to  reduce  the  distance 
from  the  surface  to  the  end  of  the  mixing  region,  measured  perpendicular  to  the 
surface,  by  che  factor  sintf,  where  0  is  the  angle  between  the  surface  and  the 
resultant  vector  of  the  transpiration  and  crossflow  velocities. 

Heat  release  from  the  AP  monopropellant  flame  is  assumed  to  occur  at  one 
plane,  resulting  in  a  discontinuity  in  the  temperature  derivative  at  its  point 
of  release,  while  the  columnar  diffusion  flame  is  assumed  to  release  its  heat 
in  a  distributed  fashion  (the  distribution  being  defined  by  a  Burke- 
Schumann^)  analysis)  between  distances  Lpy  and  Lpx  +  FH9Osin0  from  the 
surface. 

Details  of  the  equation  development  for  the  unimodal  oxidizer,  non- 
metalized  propellant  model  and  the  solution  procedure  are  presented  in 
References  1  and  2.  Included  in  the  model  are  three  "free"  constants  to  be 
chosen  by  optimization  against  data.  These  three  constants  are  pre¬ 
exponentials  associated  with  the  subsurface  heat  release  rate  expression  and 
the  two  rate  expressions  used  to  calculate  the  gas-phase  reaction  times.  These 
constants  were  chosen  on  the  basis  of  no-crossflow  burning  rate  versus  pressure 
data  for  four  unimodal  oxidizer  AP/HTPB  (hydroxy terminated  polybutadiene) 
formulations  and  then  used  unchanged  for  no-crossflow  and  crossflow  calcu¬ 
lation  for  multimodal  oxidizer  and  metalized  formulations. 
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Extension  of  the  basic  model  to  treat  multimodal  oxidizer  propellants  was 
carried  out  in  a  very  straightforward  manner  using  a  minor  variation  on  Glick's 
"petit  ensemble"  approach^’).  In  Glick's  approach,  a  propellant  containing 
oxidizer  particles  of  different  sizes  is  broken  into  a  series  of  sub- 
prop'el lants ,  each  of  which  contains  oxidizer  of  only  .  1  size.  The  sub¬ 
propellents  are  assumed  to  burn  non-interactively ,  with  tiie  unimodal  oxidizer 
model  being  used  to  calculate  a  mass  flux  for  each,  and  straightforward 
averaging  weighted  according  to  fractional  surface  areas  associated  with  each 
subpropellant  then  being  used  to  obtain  an  overall  propellant  average  linear 
regression  rate.  The  only  manner  in  which  oxidizer  of  one  size  is  allowed  to 
affect  the  burning  rate  of  a  subpropellant  containing  oxidizer  of  another  size 
is  through  possible  influence  on  the  assignment  of  fuel  to  that  subpropellant. 
That  is,  rather  than  fuel  having  to  be  assigned  to  each  oxidizer  size  category 
in-direct  proportion  to  the  amount  of  oxidizer  in  that  category,  the  capability 
of  allowing  uneven  assignment  of  fuel  to  various  oxidizer  size  subpropellants 
is  included  by  means  of  a  power  latv : 

Vf.di  ■  C2  (D0)xexp°F  (1) 

where  Vj  d-  is  the  volume  of  fuel  assigned  to  a  particle  of  diameter  (D0)f, 
XEXPOF  is  an  arbitrary  input  power  law  constant,  and  C2  is  a  constant 
determined  by  application  of  overall  continuity.  It  may  easily  be  shown  that 
XEXPOF  =  3  will  result  in  each  subpropellant  having  the  same  oxidizer/fuel 
ratio  as  the  overall  propellant.  XEXP0F<3,  on  the  other  hand,  will  result  in 
subpropellants  with  small  oxidizer  being  more  fuel-rich  than  the  overall 
propellant  and  subpropellants  with  large  oxidizer  being  more  fuel-lean,  with 
the  reverse  occurring  for  XEXP0F>3.  (Based  on  preliminary  calculations,  XEXPOF 
has  been  set  equal  to  3  in  the  calculations  presented  herein.) 

As  mentioned,  one  modification  to  the  Glick  approach  has  been 
introduced.  His  averaging  procedure  involves  an  implicit  assumption  that  the 
average  fractional  surface  area  of  the  total  propellant  which  is  subpropellant 
i  is  equal  to  the  volume  fraction  of  the  propellant  which  is  subpropellant  i. 
However,  careful  examination  of  the  situation  indicates  that  if  there  are 
subpropellants  burning  at  different  rates,  the  slower  burning  ones  will  at  any 
given  time  occupy  a  higher  fraction  of  the  surface  than  indicated  by  their 
volume  fraction:  thus,  it  seems  that  an  averaging  approach  based  on  residence 
time  distributions  is  more  appropriate.  Development  of  this  concept  leads  to 
the  average  burning  rate  for  XEXP0F=3  cases  being  given  by: 

r  =  1/Sxi/ri  (2) 

rather  than  by  Glick's  formula: 

r=vxiri  •  (3) 

ri  =  burning  rate  of  subpropellant  i 
xi  =  mass  fraction  of  propellant  in  subpropellant  i 
r  =  average  burning  rate 

The  model  described  in  detail  in  Reference  2  has  recently  been  extended  to 
treat  the  effects  of  aluminum  additive  on  propellant  burning  rate.  The 
metal  is  allowed  to  affect  burning  rate  both  through  heat  sink  effects  (altered 
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heatup  terms  in  the  overall  energy  balance  and  a  depressed  gas  diffusion  flame 
temperature  relative  to  that  which  would  be  achieved  by  the  oxidizer  and  fuel 
in  the  absence  of  the  metal  heat  sink)  and  via  conductive  and  radiative 
feedback  of  heat  released  by  the  metal  particles  burning  with  the  gas  flame 
products  above  the  propellant  surface.  Among  the  phemonema  treated  in  this 
model  are  aluminum  agglomeration  at  the  surface,  particle  velocity  lag 
relative  to  the  gases  leaving  the  propellant  surface,  particle  ignition  delay, 
particle  combustion,  conductive  feedback  from  incremental  heat  release  zones 
at  various  distances  from  the  propellant  surface,  and  radiative  feedback.  In 
the  case  of  multimodal  oxidizer  propellants,  the  assignment  of  various 
fractions  of  the  aluminum  to  the  various  subpropellants  is  treated  in  the  same 
manner  as  the  assignment  of  fuel  (binder)  to  these  subpropellants.  (As  with 
the  assignment  of  binder  to  subpropellants  containing  different  size  oxidizer 
particles,  in  the  calculations  presented  herein,  the  assignment  of  metal  was 
carried  out  such  that  all  subpropellants  were  of  the  same  composition  as  the 
overall  formulation.) 


Aluminum  particles  are  assumed  to  leave  the  propellant  surface  in  two 
forms,  agglomerated  and  unagglomeratod  (virgin),  with  the  agglomerate  fraction 
and  size  being  calculated  using  empirical  formulae  developed  by  Beckstead^ . 


Dag  =  max  (°un>  (i-()  —  1  ly,3Dox,c) 

Q  nV  r 1  * 


f  -  1-0.0005  Dqx> c  dun 


0.25 


ox,c 


I t+30(F/C) 
[l+80(F/C)2 


(4) 

(5) 


UAG 

dUN 

D0X,C 

^A1 

<’OX,C 

F/C 


=  Fraction  aluminum  not  agglomerated 

r  Agglomerate  diameter  (microns) 

=  Virgin  metal  diameter  (microns) 

=  Coarse  fraction  oxidizer  diameter  (microns) 

=  Volume  fraction  Aluminum  in  Propellant 

=  Volume  Fraction  Coarse  Oxidizer  in 
Propellant 

=  Weight  of  Finest  Oxidizer  Cut/Weight  of 
Coarsest  Oxidizer  Cut 


Based  on  preliminary  calculations,  agglomerates  were  assumed  to  con¬ 
tribute  negligible  conductive  feedback  to  the  proepllant  surface  due  to  long 
ignition  and  combustion  times.  Ignition  delay  times  for  virgin  particles 
subsequent  to  leaving  the  surface  are  calculated  from  a  simple  heatup  formula: 


t .  =  PA1 

tlsn  TTa 


JIM 


gns 


cp,Al  *n 


/t Final~Ts  \+  hue  I  t 
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=  Ignition  delay  time 
=  Aluminum  particle  density 
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Agas  =  8as  thermal  conductivity 

Cp,Al  =  Aluminum  specific  heat 

Tp£nai  =  Equilibrium  flame  temperature  subsequent  to  A1 
combustion 

Ts  =  Propellant  Surface  Temperature 

Tj&;<  =  Aluminum  particle  ignition  temperature  (2100°K) 

^melt  =  Aluminum  melting  temperature 

^melt  =  Aluminum  heat  of  fusion 

After  ignition,  the  particle  diameter-time  history  (and  thus  burning 
rate)  is  calculated  using  a  rearranged  form  of  the  Belyaev  (8)  burning  time 
formula : 


Ot+At  = 

[Dt1-5  -  (ak°-9/670)  At j  °-667 

Dt 

-  particle  diameter  (cm)  at  time  t 

Dt+At 

-  particle  diameter  at  time  t+At 

At 

=  time  increment  (seconds) 

ak 

=  100  (nC02  +  nH?o  +  n02)/-ni 

Thus,  from  Equations  6  and  7,  the  fractional  consumption  of  aluminum  as  a 
function  of  time  subsequent  to  its  leaving  the  surface  is  calculated.  (Note 
that  implicit  in  Equation  6  is  an  assumption  that  the  aluminum  particle 
temperature  is  equal  to  the  propellant  surface  temperature  when  it  leaves  the 
surface . ) 

For  calculation  of  the  distance  of  the  particle  from  the  surface  versus 
time  (required  for  calculation  of  the  metal  heat  release  distribution)  the  gas 
velocity  component  normal  to  the  surface  is  first  calculated  as  a  function  of 
distance  from  the  surface  using  the  propellant  burning  flux  calculated  in  the 
previous  loop  of  the  overall  trial  and  error  solution  procedure  along  with 
assumption  of  a  linear  temperature  profile  from  the  propellant  surface  to  the 
end  of  the  columnar  diffusion  flame  zone  with  subsequent  constant  temperature. 
Next,  a  force  balance  on  the  particle  (with  a  buoyancy  term  neglected  on  the 
basis  of  preliminary  estimates  for  particle  sizes  of  interest): 

Pg  CD  (vg~vp)2  (8) 

/>A1  D 

=  change  in  particle  velocity  in  time  increment  At 
=  gas  density 
=  drag  coefficient 
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Vg  =  gas  velocity  component  normal  to  surface 

Vp  =  particle  velocity  component  normal  to 

surface 

D  =  particle  diameter 

is  integrated  outward  from  the  surface,  utilizing  a  fitted  drag  coefficient 
versus  Reynolds'  Number  expression,  with  this  coefficient  assumed  (based  on 
the  work  of  Marshall  (9))  to  increase  by  a  factor  of  2.5  for  ignited  particles. 
(Actually,  Equations  7  and  8  have  to  be  integrated  simultaneously  subsequent  to 
particle  ignition  since  they  are  coupled  through  the  changing  particle 
diameter.)  With  coupling  of  Equations  5-8,  the  fraction  of  the  total  aluminum 
in  the  propellant  burned  is  calculated  as  a  function  of  distance  from  the 
propellant  surface.  This  is  then  multiplied  by  the  heat  of.  combustion  of 
aluminum  in  the  oxidizer/fuel  products,  calculated  as: 

QAL  -  cpm  (TFinal-T*)/XAl  '  (9) 

QaI  =  aluminum  heat  of  combustion  per  unit  mass 

\ 

Cpm  =  mixture  specific  heat 

T*  =  equilibrium  flame  temperature  with 

non-reacting  aluminum 

XAi  =  weight  fraction  aluminum  in  the 
propellant 

to  yield  a  distribution  of  heat  release  from  aluminum  combustion  versus 
distance  from  the  propellant  surface.  Conductive  feedback  flux  from  the 
aluminum  combustion  to  the  propellant  surface  is  then  calculated  from: 

qfdbk.Al.c  -  Qj  TT  expl-mcpAxi/Ai)  <10) 

j  *  i  *1 

m  =  propellant  burning  mass  flux  (gm/cm^sec) 

qfdbk,Al,c  =  aluminum  feedback  flux  cal/cm^sec) 

Qj  =  aluminum  heat  release  per  unit  mass  of 

propellant  at  jth  node  from  the 
surface  (cal/gm) 

Cp  =  product  heat  capacity 

=  size  of  ith  increment 

=  total  thermal  conductivity  (laminar 
plus  turbulent) associated  with 
ith  increment 


UNCLASSIFIED 


UNCLASSIFIED 


Preliminary  calculations  indicated  that  radiative  feedback  should  in 
general  be  small,  but  a  very  simplified  analysis  for  estimating  this  feedback 
was  included.  In  this  analysis,  an  effective  distance  from  the  surface  is 
input  as  a  free  parameter  and  an  effective  cloud  emmisivity  based  on  the 
fractional  area  subtended  by  particles  in  that  region  is  then  calculated:  for 
motors,  the  distance  is  chosen  to  be  the  port  radius  while  for  strands  it  is 
arbitrarily  set  equal  to  1  to  2  centimeters  (with  the  assumption  that  cooling 
due  to  entrained  nitrogen  from  the  surroundings  makes  particles  further  away 
ineffective).  For  the  metalized  formulations  tested  and  discussed  in  this 
paper,  radiative  transfer  is  predicted  to  have  an  insignificant  effect  on 
predicted  burning  rate. 

It  was  originally  postulated  by  this  author  that  the  augmentation  of 
composite  propellant  burning  rate  due  to  crossflow  resulted  solely  from 
shortening  of  the  distance,  measured  normal  to  the  propellant  surface, 
associated  with  mixing  of  the  fuel  and  oxidizer  gases  through  distortion  of  the 
mixing  cone,  as  depicted  in  Figure  2.  A  detailed  discussion  of  this  mechanism 
is  given  in  Reference  10.  From  tail-flame  theory,  for  no  crossflow,  the  mixing 
distance  is  related  to  the  burning  mass  flux  and  oxidizer  particle  diameter  by: 

FH90  =  k  tiipdp2  (11) 

With  crossflow  (refer  to  Figure  2)  the  mass  flux  in  the  direction  of  the 
resultant  flow  is: 

m  =  mp/sintf  (12) 

where  6  is  the  resultant  flow  angle.  However,  the  characteristic  mixing  time 
is  decreased  since  the  average  mixing  concentration  gradient  is  increased  by 
the  circular  mixing  cross-section  (in  the  absence  of  crossflow)  being 
converted  to  an  elliptical  cross-section  with  major  axis  dp  and  minor  axis 
dpsinG.  An  exact  mixing  calculation  for  this  geometry  isHof  course  quite 
difficult,  but  replacement  of  the  circle  diameter  dp  by  the  geometric  mean 
ellipse  diameter,  Vdp-sin0  in  Equation  11  does  not  seem  unreasonable.  With 
substitution  of  (12)  into  (11),  one  finds: 

krn 

FH90  ~  — (d_2  sin0)  =  kmD  dD2  (13) 

Sin a  r  r  r 

as  in  the  case  of  no  crossflow.  However,  this  FH90  is  now  measured  along  the 
resultant  flow  vector:  accordingly,  the  distance  from  the  surface  to  the  end  of 
the  mixing  zone  is  reduced  to  FH9Osin0. 

Use  of  an  erosive  burning  package  based  solely  on  this  mechanism  in  the 
different  Generation  2  model  variants  led  in  all  cases  to  major  under- 
prediction  of  the  effect  of  crossflow  on  burning  rate,  indicating  that  the 
flame-bending  mechanism  was  by  itself  insufficient.  Accordingly,  a  second 
possible  mechanism,  augmentation  of  turbulence  transport  properties  in  the 
region  between  the  propellant  surface  and  the  gas-phase  flames  was  invoked  and 
combined  with  the  flame-bending  mechanism.  In  this  approach,  it  was  assumed 
that  both  the  effective  thermal  conductivity  (governing  feedback  from  the 
various  gas  flames)  and  the  effective  mass  diffusivity  (an  important  parameter 
in  determining  the  thickness  of  the  diffusion  flame)  were  increased  in 
crossflow  situations  by  crossflow-induced  turbulence.  A  flow  profile  analysis 
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permitting  calculation  of  eddy  viscosity  (and,  by  analogy,  total  effective 
thermal  conductivity  and  diffusivity)  as  a  function  of  distance  from  the 
propellant  surface  for  a  given  crossflow  velocity,  transpiration  velocity 
(determined  by  the  propellant  burning  rate),  and  temperature  field  (depending 
on  the  location  of  gas-phase  heat  release  zones)  was  developed  and  coupled  with 
the  combustion  model  for  erosive  burning  calculations.  An  improved  cal¬ 
culation  of  diffusion  flame-bending  angle  was  also  incorporated  in  this 
analysis. 


Details  of  the  flow  profile  analysis  procedures  are  presented  in 
Reference  2.  (A  summary  of  the  equations  used  and  the  solution  procedure 
employed  are  given  in  Table  I  and  II.)  The  outputs  from  this  analysis  are  used 
to  calculate : 


^ef fective/^laminar  =  ^ef fective/^laminar  -  1  +pt/ti~  F(y) 


That  is,  the  ratio  of  total  transport  properties  to  laminar  transport 
properties  are  calculated  as  a  function  of  distance  from  the  surface.  Average 
total  transport  property  values  between  appropriate  zones  are  then  calculated 
and  substituted  for  the  laminar  values  in  the  diffusional  mixing  equations  and 
the  heat  feedback  equations  in  the  original  model,  revised  burn  rates  and  flame 
distances  are  calculated,  and  the  procedure  is  repeated  until  convergence  is 
achieved.  As  might  be  expected,  this  looping  procedure  is  considerably  more 
complex  in  the  case  of  multimodal  propellants  than  for  unimodal  propellants 
since  solution  of  the  individual  subpropellant  cases  becomes  interactive  in 
the  case  of  crossflow.  This  interaction  occurs  because  there  is  only  one 
boundary  layer  for  the  overall  propellant-  (that  is,  one  cannot  calculate  a 
different  boundary  layer  profile  for  each  subpropellant)  with  the  boundary 
layer  details  being  controlled  by  the  average  transpiration  velocity,  flame 
height,  surface  temperature,  etc.  for  the  overall  propellant  rather  than  by  the 
individual  values  of  these  parameters  for  each  subpropellant. 


It  should  be  noted  that  the  procedures  used  for  calculating  average 
transport  properties  between  zones  from  the  distribution  of  these  properties 
have  been  modified  from  those  procedures  used  in  calculation  of  predicted 
burning  rates  presented  in  References  2  and  11.  In  the  earlier  work,  linear 
averaging  was  employed,  while  the  modified  procedure  allows  for  the  fact  that 
the  effect  of  the  gas-phase  heat  release  on  burning  rate  decreases  exponen¬ 
tially  with  the  distance  of  that  release  from  the  surface.  In  addition,  in 
calculation  of  results  presented  previously,  the  bent  flame  was  approximated 
as  being  straight  with  the  effective  angle  being  determined  by  integration  of 
the  velocity  vector  profile  from  the  surface  to  the  end  of  the  diffusion  flame 
zone  while  the  effects  of  flame  curvature  are  included  in  the  current  model. 


Several  options  for  closure  of  the  boundary  layer  analysis  through 
specification  of  an  eddy  viscosity  equation  were  built  into  the  model:  all 
entailed  use  of  a  Prandtl  mixing-length  expression  of  the  general  form: 


<  =  .168  (y+yA)2  (DF)2  du/dy 


where  DF  is  a  damping  factor  which  is  a  function  of  such  parameters  as  blowing 
ratio,  axial  pressure  gradient,  and  roughness  height  while  yA  is  an  offset 
factor  dependent  on  roughness  height.  Most  satisfactory  results  were  found 
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with  DF  set  equal  to  unity  (no  damping)  and  set  equal  to  zero  (no  roughness 
effect). 


EXPERIMENTAL 

The  experimental  fest  apparatus  (Figure  3)  and  procedures  for  measuring 
burning  rates  with  crossflow  are  described  in  detail  in  Reference  12.  A 
cylindrical ly  perforated  "driver"  grain  whose  length  is  chosen  eo  give  the 
desired  operating  pressure  produces  a  high  velocity  gas  flow  through  a 
transition  section  into  a  rectangular  thannel  which  contains  the  test  grain. 
The  test  grain  extends  through  the  transition  section  to  butt  against  the 
driver  grain  to  eliminate  leading  edge  effects.  The  test  grain  is  ap¬ 
proximately  40  cm  long,  1.90  cm  wide,  and  2.50  cm  deep  (web)  and  burns  only  on 
the  1.90  cm  face.  The  flow  channel  is  initially  1.90  x  1.90  cm,  opening  up  to 
1.90  cm  x  4.40  cm  as  the  test  propellant  burns.  For  high  crossflow  velocity 
tests,  the  apparatus  is  operated  without  a  nozzle  while  for  lower  velocity 
tests,  a  two-dimensional  nozzle  is  employed. 

The  burning  rate  is  directly  measured  by  photographing  the  ablating  grain 
with  a  high-speed  motion  picture  camera  through  a  series  of  four  quartz  windows 
located  along  the  length  of  the  test  section.  Frame  by  frame  analysis  of  the 
films  permits  determination  of  instantaneous  burning  rate  as  a  function  of  time 
at  each  of  the  four  window  locations. 

For  nozzled  cases,  the  measured  location  of  the  burning  propellant 
surface  at. each  window  as  a  function  of  time,  together  with  the  known  constant 
throat  area,  permits  straightforward  calculatio  .  of  the  crossflow  velocity  as 
a  function  of  time.  However,  the  very  sensitive  dependence  of  Mach  number  on 
area  ratio  for  M>0.5  makes  calculation  of  crossflow  velocity  from  area  ratio 
measurement  quite  poor  for  nozzleless  cases.  Accordingly,  for  these  tests, 
stagnation  pressure  is  determined  at  the  aft  end  of  the  test  section  and  used 
in  combination  with  the  driver  chamber  pressure  for  calculation  of  the 
stagnation  pressure  in  the  test  section  as  a  function  of  time  and  position. 
Static  presurc  wall  taps  at  each  window  location  are  used  for  measurement  of 
static  pressure  as  a  function  of  time  for  both  nozzled  and  nozzleless  cases. 
From  the  static  and  stagnation  pressure  values  determined  as  a  function  of  time 
and  position  down  the  test  section,  crossflow  Mach  Number  and  velocity  are 
calculated  as  a  function  of  time  at  each  window  location  in  the  test  section  for 
the  nozzleless  cases. 

The  erosive  burning  characteristics  of  a  series  of  11  AP/HTPB  formula¬ 
tions  (4  with  unimodal  AP  and  no  catalyst  or  metal,  5  with  multimodal  AP  and  no 
catalyst  or  metal,  1  with  multimodal  AP  and  aluminum  but  no  catalyst,  and  1  with 
unimodal  AP  plus  catalyst)  with  systematically  varied  parameters  have  been 
measured  in  this  device  to  date.  In  addition,  standard  strand-burning 
procedures -have  been  used  to  determine  zero-crossflow  burning  rate  versus 
pressure  character ist ics  for  these  11  formulations  plus  3  additional  multi¬ 
modal  AP,  aluminized  (no  catalyst)  formulations.  The  propellant  matrix  tested, 
including  the  rationale  for  its  selection,  is  presented  in  Table  III. 

EXPERIMENTAL  RESULTS  AND  COMPARISON  WITH  MODEL  PREDICTIONS 

Predicted  and  measured  zero-crossflow  burning  rate  versus  pressure 
characteristics  of  13  of  the  14  formulations  are  presented  in  Figure  4-6. 
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(Predictions  were  not  made  Cor  the  catalyzed  formulation  since  the  catalyst 
should  cause  shifting  of  one  or  more  of  the  three  aforementioned  kinetic 
constants  and  it  was  felt  that  the  catalyzed  formulation  data  base  was 
insufficient  for  re-evaluation  of  these  constants.)  As  may  be  seen,  agreement 
between  predicted  and  measured  200  psia  burning  rates  is  excellent,  with  all 
predictions,  covering  a  wide  range  of ‘ non-metalized  and  metalized  AP/HTPB 
formulations,  falling  within  10  percent  of  the  measured  values.  Similarly,  the 
1000  psia  burning  rate  predictions  and  data  show  excellnt  agreement,  with  only 
one  point  falling  outside  the  20  percent  error  band  and  one  falling  between 
the  10  and  20  percent  bands,  the  remainder  agreeing  within  10  percent.  As 
indicated  by  Figure  6,  the  1000  psia  predicted  versus  measured  value 
comparisons  of  burning  rate  pressure  exponent  are  also  quite  good,  all  lying 
within  the  20  percent  error  bands  and  most  within  the  10  percent  bands. 
Accordingly,  it  appears  that  development  of  an  accurate  zero-crossflow  model , 
necessary  as  a  first  step  in  developing  an  adequate  erosive  burning  model,  has 
been  satisfied. 

Erosive  burning  test  results  are  presented  in  Figures  7-17,  in  the  form  of 
burning  rate  versus  pressure  at  various  crossflow  velocities.  In  addition, 
theoretical  predictions  are  presented  for  all  but  the  catalyzed  formulation 
(4869)  in  these  figures.  As  may  be  seen,  the  agreement  between  data  and 
predictions  for  the  no-crossflow  conditions  is  excellent,  as  indicated 
earlier,  for  all  formulations  except  7996,  where  the  theoretically  predicted 
burning  rates  are  10  to  20  percent  high.  In  addition,  the  crossflow  effect 
predictions  agree  reasonably  well  with  the  data  in  general.  With  the  baseline 
formulation  (4525),  the  theory  slightly  underprcdicts  the  effect  of  crossflow 
on  burning  rate  while  with  5051,  4685;  and  5542  (the  other  three  non-catalyzed 
unimodal  ozidizer  formulations)  agreement  between  theory  and  data  is  excel¬ 
lent.  The  model  predicts  that  the  high  burning  rate  formulation  (5555)  should 
be  quite  insensitive  to  crossflow  velocity,  in  excellent  agreement  with 
experiment.  Formulation  5565,  on  the  other  hand,  appears  to  be  slightly  less 
sensitive  to  crossflow  than  predicted,  particularly  at  the  lower  pressures  (1- 
3  MPa).  Agreement  between  theory  and  data  for  the  remaining  three  multimodal 
oxidizer,  non-metalized  formulations  is  good,  except  that  the  zero-crossflow 
offset  between  theory  and  data  for  7996  appears  to  be  maintained  for  the 
crossflow  cases.  Finally,  the  rather  limited  data  for  the  metalized 
formulation  (6626)  appear  to  be  in  general  agreement  with  predictions. 

Results  for  the  various  formulations  may  be  compared  to  identify 
parameters  dominating  the  sensitivity  of  burning  rate  to  crossflow.  Formula¬ 
tions  4525,  5051,  and  4685  were  identical  except  for  oxidizer  particle  size 
(and,  as  a  consequence)  base  (no-crossflow)  burning  rate.  Examination  of 
Figures  7-9  reveals  that  the  crossflow  sensitivity  increases  with  increasing 
particle  size  (decreasing  base  burning  rate).  '  For  example,  at  200  m/scc  and  5 
MPa,  the  augmentation  ratios  for  4685,  4525  and  5051  are  about  1.10,  1.60  and 
2.00  respectively. 

Comparison  of  data  for  4525  and  4869,  differing  only  in  use  of  catalyst  in 
the  latter  (with  consequent  higher  base  burn  rate)  again  shows  an  increase  in 
crossflow  sensitivity  with  decreasing  base  rate.  At  5  MPa  and  200  m/sec,  their 
respective  burn  rate  augmentation  ratios  are  1  .60  and  1.10,  while  at  600  m/sec, 
the  r/r0  values  are  2.3  and  1.7.  Thus,  base  burn  rate  is  seen  to  affect  erosion 
sensitivity  even  at  constant  oxidizer  size. 
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Formulations  4685  and  4869  have  approximately  the  same  base  burning  rate 
at  8  MPa  although  their  oxidizer  sizes  are  different.  Data  comparison 
indicates  that  these  formulations  have  nearly  the  same  sensitivity  to  low 
crossflow  velocities  at  8  MPa,  with  the  catalyzed  propellant  being  only 
slightly  more  sensitive  at  higher  velocities.  Thus,  it  appears  that  it  the 
base  burning  rate  rather  than  the  oxidizer  particle  size  which  dominates  the 
sensitivity  of  this  series  of  four  73/27  AP/I1TPB  formulations  to  crossflow, 
though  oxidizer  size  itself  does  appear  to  have  a  slight  additional  effect, 
crossflow  sensitivity  decreasing  with  decreasing  size  at  constant  base  rate. 

Formulation  5542  differs  from  4525  in  oxidizer/fuel  ratio  and  conse¬ 
quently  flame  temperature..  Since  oxidizer  particle  size  was  held  constant,  the 
higher  0/F  ratio  results  in  higher  base  rate  for  5542.  The  data  (Figure  7  and 
11)  indicate  t?hat  the  crossflow  sensitivity  of  5542  is  considerably  lower 
over  the  entire  range  of  conditions  studied.  Comparison  of  results  for  5565 
and  4525,  which  differ  in  0/F  ratio,  but  have  the  same  base  burning  behavior 
(due  to  compensating  AP  particle  size  differences)  indicates  that  the 
sensitivity  of  these  two  formulations  to  crossflow  is  nearly  identical. 
Accordingly,  it  may  be  concluded  that  0/F  ratio  (and  consequently  flame 
temperature)  changes  do  not  directly  effect  the  erosion  sensitivity  of  these 
forumulations,  but  only  affect  it  through  their  effect  on  base  burning  rate. 

As  indicated  in  Table  III,  Formulation  6626  (metalized)  has  nearly  the 
same  base  burning  characteristics  as  4525  and  5565  and  approximately  the  same 
flame  temperature  as  5565.  The  data  of  Figures  7,  12  and  17  reveal  that  all 
three  formulations  have  quite  similar  erosive  burning  characteristics.  For 
example,  at  760  m/sec  and  2.8  MPa,  the  augmentation  ratios  for  4525,  5565,  and 
6626  are  2.1,  2.25,  and  2.1,  while  at  260  m/sec  and  4.0  MPa,  they  are  1 .65,  1.55 
and  1.60.  These  results  support  a  conclusion  that  the  dominant  factor 
affecting  crossflow  sensitivity  of  composite  propellants  is  base  burning  rate, 
largely  independent  of  the  factors  determining  that  base  rate. 

Formulations  5555,  5565,  7993,  7996,  and  8019  are  identical  in  composi¬ 
tion  (82/18  AP/HTPB)  differing  only  in  oxidizer  particle  size  blends,  which 
were  adjusted  to  give  a  range  of  base  (zero-crossflow)  burning  rate  versus 
pressure  characteristics.  In  Figure  18,  data  extracted  from  Figures  12-16  are 
plotted  in  the  form  of  burning  rate  augmentation  factor  (r/rQ)  versus  base 
burning  rate  for  three  combinations  of  pressure  and  crossflow  velocity.  As  may 
be  seen,  the  augmentation  factor  decreases  monotonically  and  fairly  smoothly 
with  increasing  base  burning  rate,  again  indicating  the  importance  of  that 
parameter  on  crossflow  sensitivity. 

A  summary  of  the  above  comparisons,  delineating  the  effects  of  various 
parameters  on  crossflow  sensitivity  of  burning  rate  is  presented  in'  Table  IV. 

MOTOR  SCALING  EFFECTS 

It  has  been  observed  in  the  past  that  erosive  burning  effects  tend  to 
diminish  with  increasing  motor  size,  all  other  parameters  being  held  constant. 
It  is  therefore  of  interest  to  ascertain  whether  the  model  described  above  can 
predict  such  a  trend.  Accordingly,  a  series  of  calculations  have  been  run  for 
equal  crossflow  velocities  and  pressures  with  different  motor  diameters. 
Typical  results  are  plotted  in  Figures  19  and  20,  comparison  being  made  between 
motors  with  port  diameters  of  1.0  and  10  inches.  As  may  be  seen,  erosive 
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burning  o f fects  are  indeed  predicted  to  decrease  significantly  with  increasing 
motor  size.  Another  point  of  interest  may  be  made  from  Figure  20.  One  of  the 
weakest  links  in  the  use  of  the  flow  profile  analysis  routine  outlined  earlier 
for  calculation  of  the  flame  bending  angle  and  turbulent  transport  property 
distribution  is  specification  of  the  ratio  of  skin  friction  coefficient  with 
blowing  to  that  without  blowing  as  a.  function  of  the  blowing  ratio,  as 
discussed  in  detail  in  Reference  2.  Accordingly,  several  different  empirical 
expressions  have  been  examined  (See  Figure  21).  As  shown  in  Figure  20,  the 
predicted  burning  rate  augmentation  ratio  is  in  general  fairly  insensitive  to 
the  choice  of  expression  among  those  considered  as  reasonable  possibilities. 

CONCLUSIONS 

A  fundamental  composite  propellant  combustion  model  capable  of  predicting 
burning  rate  as  a  function  of  pressure  and  crossflow  velocity,  given  only 
propellant  composition  and  ingredient  particle  size  has  been  developed. 
Testing  against  zero-crossflow  data  for  a  series  of  13  non-catalyzed  AP/HTPB 
formulations  (4  containing  aluminum  additive)  indicates  that  the  model  does  an 
excellent  job  in  predicting  the  non-erosive  burning  rate  versus  pressure 
characteristics  of  sudh  formulations.  In  addition,  erosive  burning  data  have 
been  obtained  for  10  of  these  formulations:  the  model  also  does  a  good  job  in 
predicting  the  observed  erosive  burning  behavior  of  these  propellants. 

The  data  obtained  to  date  indicate  that  the  base  (no-crossflow)  burning 
rate  characteristics  of  the  propellant  have  a  predominant  affect  on  its 
sensitivity  to  crossflow,  high  burning  rate  formulations  being  considerably 
less  susceptible  to  erosive  burning  than  low  burning  rate  formulations, 
whether  the  base  burning  rate  alterations  are  produced  by  oxidizer  particle 
size  variation,  oxidizer/fuel  ratio  variation,  addition  of  metal,  or  use  of 
catalysts.  Thus  propellants  with  widely  differing  oxidizer  size  distribution, 
0/F  ratios,  metal  loadings,  etc  tend  to  show  identical  erosive  burning  behavior 
as  long  as  they  have  identical  base  (no  crossflow)  burning  rate  charac¬ 
teristics.  Oxidizer  particle  size  does  appear  to  have  some  residual  effect 
(but  only  a  slight  one)  beyond  its  effect  on  base  burning  rate,  erosion 
sensitivity  increasing  with  increasing  particle  size.  An  important  pre¬ 
liminary  conclusion  is  that  aluminum  (at  least  at  low  levels)  does  not  affect 
erosion  sensitivity  other  than  through  its  effect  on  base  burning  rate. 

The  model  also  predicts  a  decrease  in  erosive  effects  with  increasing 
motor  port  size  (at  equal  pressure  and  crossflow  velocity)  in  agreement  with 
observations  by  other  investigators. 
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Table  I.  Main  Equations  used  in  Flow  Profile  Analysis  (Coupled  with  Burning  Rate 
Analysis  Through  Temperature  Dependence  of  Density  and  Viscosity). 

.1.  Cf0  =  f-j  (Re,De/k) 

2.  Cf/Cf0  =  f2  (B)  B  ='  (mjnj/mcrossf,ow)/(C|0/2) 

3.  rwa||  -  (C(/2)  pfs  U* 

4.  Ty  -  rwa||  +  ihjnj  U  -  Kmom  jnt  Y 
•  5.  ry  =  (u  +  pe)  dU/dY 

6.  €  =  [f3  (Y,  minj,  dP/dx,  k)]  dU/dY 

FOR  ZERO  ROUGHNESS  AND  NO  DAMPING.  FOR  EXAMPLE 
e  =  .168  Y2  (dU/dY) 

7.  p  =  P{MW)/RT 

8.  n  =  kT0  8 

9.  T  =  Ts  +  (Tf  -  Ts)Y/(LRx  +  FH90sinO) 

10.  V  =  mjnj/p 


Table  II.  Application  of  Flow  Profile  Analysis  Equations. 

1.  USE  EONS  1-4  TO  OBTAIN  t  =  g(u,y) 

2.  ’INTEGRATE  EQN  5.  WITH  USE  OF  EONS  4  AND  6-9,  FROM  THE  WALL 
(SURFACE)  TO  OBTAIN  u.  du/dy,  e  AS  FUNCTION  OF  y 

3.  KNOWING  c(y),  COMBINE  THISWITH  EONS  7-9  TO  OBTAIN  (1  * pd») 

AS  A  FUNCTION  OF  y 

4.  USE  THIS  FUNCTION  TO  OBTAIN  AVERAGE  AUGMENTATION  RATIOS 
FOR  LAMINAR  TRANSPORT  PROPERTIES  OVER  VARIOUS  REGIONS  AND 
USE  THESE  TO  ADJUST  EFFECTIVE  FLAME  OFFSET  DISTANCES 

5.  HAVING  u(y)  AND  w(y)  INTEGRATE  ALONG  A  STREAMLINE  TO  DETERMINE 
THE  EFFECTIVE  REDUCTION  OF  DIFFUSIONAL  DISTANCE  DUE  TO  FLAME 
BENDING. 

6.  LOOP  -  FAIRLY  COMPLEX  FOR  UNIMODAL  AP  CASES;  EVEN  MORE  SO  FOR 

'  MULTIMODAL  AP  CASES.  WHERE  INTERACTION  BETWEEN  THE  MODES  MUST 
NOW  BE  CONSIDERED 


•FOR  LOW  CROSS  FLOW.  HIGH  BLOWING  RATIO  CASES,  THE-Ky  TERM  IN  EQN4 
CAUSES  r  TO  GO  NEGATIVE,  A  RESULT  WE  INTERPRET  AS  INDICATING  BOUNDARY 
LAYER  BLOWOFF  -  IF  THIS  OCCURS,  WE  USE  A  COSINE  LAW  VELOCITY  PROFILE 
FOR  OUR  ANALYSIS  IN  PLACE  OF  INTEGRATION  OF  EQN  5. 
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Table  IV. 


Table  III.  Propellant  Matrix  (AP/HTPB)  Tested. 

FORMULATION  COMPOSITION  RATIONALE 


4525 

73/27  AP/HTPB.  20  MICRON  AP 

5051 

73/27  AP/HTPB.  200  M'CRON  AP 

4685 

73/27  AP/HTPB.  5  MICRON  AP 

4869 

72/26/2  AP/HTPB/ F*203. 

20  MICRON  AP 

5542 

77/23  AP/HfPB.  20  MICRON  AP 

•  5565  ' 

82/18  AP/HTPB.  13.65%  90  MICRON 
AP. 68.35%  200  MICRON  AP 

BASELINE  FQMULATION,  FLAME 
TEMPERATURE  •  1667eK 

COMPARE  WITH  4525  FOR 

AP  SIZE  EFFECT  AND  BASE  BURNING 

RATE  EFFECT 

COMPARE  WITH  4525  AND  5051 
FOR  APSI2E  EFFECT  AND  BASE 
BURNING  RATE  EFFECT 

COMPARE  WITH  4525  FOR  BASE 
BURNING  RATE  EFFECT  AT  CONSTANT 
AP  SIZE 

COMPARE  WITH  4525  FOR  MIXTURE 
RATIO  ANO  FLAME  TEMPERATURE  EFFECT 
AT  CONSTANT  AP  SIZE.  T*206SCK 

AP  SI2ES  CHOSEN  TO  MATCH  BASF. 
BURNING  RATE  OF  4525.  COMPARE 
WITH  4525  FOR  MIXTURE  RATIO  ANO 
FLAME  TEMPERATURE  EFFECT.  T-2575°K 


5555  82/18  AP/HTPB.  41%  1  MICRON 

AP.  41%  7  MICRON  AP 

7993  82/18  AP/HTPB.  41%  7  MICRON 

AP,  41%  90  MICRON  AP 

7996  82/18  AP/HTPB.  41%  20  MICRON 

AP.  41%  200  MICRON  AP 

8019  82/18  AP/HTPB.  27.3%  1  MICRON 

AP.  27.3%  20  MICRON  AP. 

27.4%  200  MICRON  AP 

6626  74/21/5  AP/HTPB/AI.  70% 

90  MICRON  AP.  4%  200  MICRON  AP 


8495  73/17/10  AP/HTPB/AI.  7.3% 

90  MICRON  AP.  58  4%.  200 
MICRON  AP.  7.3%.  400  MICRON 
AP 

8497  80/15/5  AP/HTPB/AI.  8  0%  90  * 

MICRON  AP.  64.0%.  200  MICRON 
AP. 8  0%  400  MICRON  AP 

WC19M  70/12/18  AP/HTPB/AI.  21.0%  20 

MICRON  AP.  39.0%.  90  MICRON 
AP.  10.0%.  200  MICRON  AP 


COMPARE  WIT*  I  5565  FOR  EFFECT 
OF  BASE  BURNING  RATE. 


\  FURTHER  STUDY  OF  APSI2E 
f  AND  BASE  BURNING  RATE  EFFECTS 


SAME  FLAME  TEMPERATURE  AND  BASE 
BURNING  RATE  AS  5565.  COMPARE 
WITH  5565  FOR  Al  EFFECT. 


STUDY  OF  EFFECTS  OF  VARIOUS 
\  ALUMINUM  LOADINGS  ANO  O/F 
f  RATIOS  ON  EROSIVE  BURNING 
BEHAVIOR. 


Effects  of  Various  Parameters  on  Sensitivity  of  Formulations  to  Crossflow. 


COMPARISON 

PARAMETERS  STUDIED 

EFFECT  ON  EROSIVE  BURNING 

4525.  5051.  46(5 

VARIED  4,.  r.  AT  FIXED  BIRDER  TYPE.  FIXEO  FLAME 
TEMPERATURE 

4, 

4524.  4(69 

VARIEO  r#  AT  FIXEO  AP  SIZE.  BIROER  TYPE.  AND 

FLAME  TEMPERATURE 

-.4  ♦£* 

4(65. 4(69 

VARIEO  dp  AT  FIXEO  8JN0CR  TYPE.  AHO  FLAME 
TEMPERATURE 

(fp  j  -*-£^suchuy 

4525. 5542 

VAR10  OF  RATIO  (AHO  THUS  FLAME  TEMPERATURE) 

AHO  r9  AT  FIXED  8IH0ER  TYPE  AHO  FIXEO  AP  SIZE 

t,4  4-*-£l 

3565. 4525 

VARIEO  O  f  RATIO  (ANO  THUS  FLAME  TEMPERATURE) 

AT  FIXEO  8W0ER  TYPE  ANO  FIXEO  r# 

T,4  ♦  £  WCHANGEO 

55(5. 3335.  2335 

2596.  Mil 

VARIEO  d„.  r.  AT  FIXEO  BINDER  TYPE.  FIXEO  FLAME 
TEMPERATURE 

35(3.(626 

ALUMINUM  VERSUS  NON  ALUMINUM  AT  FIXEO  rf 

IMOER  TYPE,  ANO  FLAME  TEMPERATURE 

£  UNCHANGE 0 

coNtiuws  moM  «aov£  cohnanisons 

1.  IHt  AUGMENTATION  TACTON  IS  SINONGlY  0EBENDTN1  ON  BASE  BURKING  NATE. 

I.  THEffi  6  A  IMAM.  AESIOUAl  EFFECT  OF  0  JUMEA  BAATCIE  SITE.  AT  HIED  IUNKINC  NATE. 

I  0.1  (UAME  TEMNENATUNCl  Ef  FECTS  £  FOB  HTNI  SYSTEMS  OXIY  THAOVGM  EFFECT  ON  IASE  AUNNMC  NATE. 
4.  *1  FIXEO  IASE  lUNNMt  NATE.  AlUMMUM  NAS  NO  EFFECT  ON  C . 
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Figure  3.  Schematic  Drawing  of  Erosive  Burning  Test  Apparatus. 
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.  Figure  4.  Zero-Crossflow  Burning  Rates  at  200  psia.  Data  and  Predictions. 
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Figure  7.  Burning  Rate  Predictions  (Solid 
Lines)  and  Data  (Points)  for 
Formulation  4525  (73/27 
AP/HTPB,  20  MICRON  AP). 
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Figure  8.  Burning  Rate  Predictions  (Solid 
Lines)  and  Data  (Points)  for 
Formulation  5051  (73/27 
AP/HTPB,  200  Micron  AP). 
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Figure  9.  Burning  Rate  Predictions  (Solid 
Lines)  and  Data  (Points)  for 
Formulation  4685  (73/27 
AP/HTPB,  5  Micron  AP). 
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Figure  10.  Burning  Rate  (No  Predictions) 
for  Formulation  4869  (72/26/2 
AP/HTPB/Fe203.  20  Micron 
AP). 
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Figure  1 1.  Burning  Rate  Predictions  (Solid 
Lines)  and  Data  (Points)  for 
Formulation  5542  (77/23 
AP/HTPB,  20  Micron  AP). 


Figure  12.  Burning  Rate  Predictions  (Solid 
Lines)  and  Data  (Points)  for 
Formulation  5565  (82/18 
AP/HTPB,  13.65%  90  Micron 
AP,  68.35%  200  Micron  AP). 
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Figure  13.  Burning  Rate  Predictions  (Solid 
Lines)  for  Formulation  5555 
(82/18  AP/HTPB,  41%  1 
Micron  AP,  41%  7  Micron  AP). 


Figure  14.  Burning  Rate  Predictions  (Solid 
Lines)  and  Data  (Points)  for. 
Formulation  7993  (82/18 
AP/HTPB,  41%  7/rAP,  41% 

90/r  AP). 
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Figure  15.  Burning  Rate  Predictions  (Solid 
Lines)  and  Data  (Points)  for 
Formulation  7996  (82/18 
AP/HTPB,  41%  20/iAP,  41% 
200/r  AP). 
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Figure  16.  Burning  Rate  Predictions  (Solid 
Lines)  for  Formulation  8019 
(82/18  AP/HTPB,  27.3% 

1/i  AP,  27.3%  20/iAP,  27.4% 
200m  AP). 


Figure  17.  Burning  Rate  Predictions  (Solid  Figure  18,  Summary  of  Results  for  82/18 
Lines)  and  Data  (Points)  for  AP/HTPB  Formulations. 

Formulation  6626  (74/21/5 
AP/HTPB/AI,  5  Micron  Al, 

70%  90  Micron  AP,  4%  200 
Micron  AP). 
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